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ABSTRACT 

The initial assembly of Space Station Freedom involves the Space Shuttle, its 
Remote Manipulation System (RMS) and the evolving Space Station Freedom. The 
dynamics of this coupled system involves both the structural and the control system 
dynamics of each of these components. The modeling and analysis of such an assembly is 
made even more formidable by kinematic and joint nonlinearities. 

The current practice of modeling such flexible structures is to use finite element 
modeling in which the mass and interior dynamics is ignored between thousands of nodes, 
for each major component The model characteristics of only tens of modes are kept out of 
thousands which are calculated. The components are then connected by approximating the 
boundary conditions and inserting the control system dynamics. 

In this paper continuum models are used instead of finite element models because of 
the improved accuracy, reduced number of model parameters, the avoidance of model order 
reduction, and the ability to represent the structural and control system dynamics in the 
same system of equations. Dynamic analysis of linear versions of the model is performed 
and compared with finite element model results. Additionally, the transfer matrix to 
continuum modeling is presented. 
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The continuum modeling approach is seen to offer a viable alternative to finite 
element modeling. The continuum approach enables increased insight for synthesis and 
integrated control/structures design. 

NOMENCLATURE 

Symbols 

A, B, C, D 
c 

EA 

EIx, Ely 
Fo 
GA 
Gly 
F 

F 0 

I 
L 
Qu 
Qs 
s 

Q 

s 

w 

a 

Superscripts 


state vector elements, coefficients of the sinusoidal and hyperbolic 
functions 

model parameter vector 
longitudinal stiffness 
bending stiffness 
constant axial force 
lateral shear 
torsional stiffness 
force distribution function 
axial, steady force 
inertia matrix 
length of beam 
deflection coefficient matrix 
angular deflection matrix 
real part of the roots 

state vector, coefficients of sinusoidal and hyperbolic mode shape 
basis functions 
real part of root 

modal frequency, imaginary part of the roots 
angular velocity vector 


T transpose 

-1 inverse 

♦ differention with respect to t 

/ differention with respect to z 
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Subscripts 


i 

c.g. 


n 

x 

y 

z 

y 


mode index 
center of gravity 
general index 
x axis 
y axis 
z axis 

torsional axis, z 



INTRODUCTION 


The initial assembly of Space Station Freedom involves the Space Shuttle, its 
Remote Manipulation System (RMS) and the evolving Space Station Freedom. The 
assembly of the Space Station Freedom is performed by positioning and connecting 22 
modules using a remote manipulation system (RMS). The dynamics of this coupled 
system involves both the structural and the control system dynamics of each of these 
components. The numerous configurations that result from this assembly process 
necessitate an efficient procedure for accurately modeling the structural and control 
dynamics. The modeling and analysis of such an assembly is made even more formidable 
by kinematic and joint nonlinearities. 

Modeling of complex flexible spacecraft is an issue which has far reaching 
consequences in controller design and the subsequent spacecraft performance. Numerous 
difficulties in controlling flexible spacecraft have been attributed to inaccuracies in modeling 
[1]. With higher controller bandwidth, modeling issues assume greater significance. 
Increased size and more demanding control specifications promise to make high 
performance control more difficult [2], Current modeling schemes for the design and 
analysis of structural and control systems have several limitations [3]. The conventional 
approach is to use elements which are void of dynamics on the interior of their boundaries. 
The computational cost and numerical inaccuracies involved in generating solutions to these 
equations impose a practical limit to the size (and consequently the accuracy) of these 
structural dynamics models. For problems with minimal control-structure interaction, the 
finite element models are adequate. High performance control systems will however 
require increased fidelity and accuracy in the models. 

Distributed parameter modeling is proposed in this work to synthesize high fidelity 
spacecraft models. The distributed parameter models provide a single set of equations for 
control and structural dynamics. The conventional finite dimensional representation of 
complex spacecraft by the finite element method suffers from the following drawback. 

Finite element models are generally too large for control work. One performs model 
reduction to reduce the model order to controller synthesis amendable dimensions. Spill- 
over of control energy into the unmodeled modes can result in instability. The proposed 
approach represents flexible structural members by partial differential equations offering 
significant advantages in modeling, parameter estimation and the integrated design of 
control/structural systems [4], [5], [6]. The present method differs from the finite element 
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method in that an individual element can represent all the modes of that “super” element and 
produce the force and moment vectors at its boundaries. These elements are then connected 
at their boundaries to form the model of the complete structure. Bishop [7] and Snowdon 
[8] have studied applications in which a limited number of such elements have been 
connected to form simple frames. The homogenization technique [9] and [10] for repetitive 
lattice trusses is particularly useful. For spacecraft control applications, it is necessary to 
connect many distributed parameter elements to represent the structural dynamics of 
complex flexible spacecraft The software programs available for continuum modeling 
include: Poeleart’s [11] DISTEL, Taylor’s [12], PDEMOD and Anderson’s [13] BUNVIS 
program. In this work, PDEMOD is used to generate some of the results. 

This paper will discuss the generation of the system of partial differential equations 
for modeling complex, flexible spacecraft. A continuum model of the assembly 
configurations of the Space Shuttle RMS - Payload will be used to study the control 
problems involved. Continuum models are shown to have distinct advantages for control 
applications. 

This paper is organized in the following manner. The formulation of the structural 
dynamics models, the transfer matrix approach to modeling and the control system 
embedding methods are presented in the next section. The numerical results for a simple 
model of the Shuttle/RMS Payload assembly are presented. The results compare the modal 
characteristics obtained using NASTRAN with the continuum results. The concluding 
remarks identify the salient features of the proposed approach and related modeling and 
analysis accomplishments to date. 


Discussion 

The formulation of the dynamics using a set of distributed parameter elements 
connected at their boundaries is key to obtaining the objectives of optimal parameter 
estimation. The types of elements to be considered are (1) rigid body with a full inertia 
matrix, and (2) dynamic, flexible beam element. The equations of motion for each of these 
elements will be considered in turn. 



Equations of Motion 


A Newtonian or inertial frame of reference is used for the motion of all beam 
elements and rigid bodies. For example, the point of attachment in the Newtonian axis of 
reference, undeflected beam is: 



Figure 1: Diagram of a Rigid Body Attached to its 
Reference Flexible Beam 

^attach,o — ^c.g ,,o ^beam^" 

For the deflected beam: 

■^attach,/ — ^attach.O ^beam^ 

— ^c.g.,0 ^beam^* ^beam^ 

The position of the body center-of-gravity due to beam deflection is: 

^c.g„r — ^attach,/ ~ ^body r 

— ^c.g„0 ■*" ^beam r ^beam M — ^body^ 

For small angular deflections 


( 1 ) 


( 2 ) 


( 3 ) 


^body ^beam + Tilda(7^ earn W ) 


( 4 ) 






Substituting, we get 


^c.g .,1 ~ ^c.g.,0 + ^beam M Tilda (7jj eam U ) r 
= ^c.g.,0 ^beam w ^{^beam u ) 


( 5 ) 


Differentiating, we get the acceleration of the body center-of-gravity: 


Tb eam M + ^(Tbeam^ ) 


( 6 ) 


Equations of motion are written for each rigid body and the forces and moments 
imparted by the beams are taken into account In each case it is necessary to account for the 
different frames of reference and joints of attachment Equations of motion for the linear 
and angular degrees of freedom for all of the bodies are assembled into a single matrix, A. 

In the time domain the equations of motion are: 


4g=X{ Forces }/'* 

^= / body" / X{ M ° mentS } 


( 7 ) 


In the frequency domain, the linear and angular equations of motion are the basis 
for each block of elements: 


' i ' 


^Linear,./ Quj+Tj RjTjQ u j 2 ^Fi } 


\ m i m i 


^angular , J ~ Quj 7} Ij ^ ^2 ^beam-»7beain-/7Vi} 


( 8 ) 


For each case in which a rigid body has more than one beam attached, a constraint 
equation is added to the system of equations. Assembly of the equations of motion and the 
constraint equations yields the system matrix from which we get the characteristic equation: 

\A(a+Jw)\=0 /m 
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Flexible Beam Equations 


The flexible beam elements exhibit lateral bending in two axes, axial deformation, 
and torsion. The governing partial differential equations have a variety of terms so that 
parameter values can select, for example, a wave or string equation, Euler beam equation, 
or Timoshenko beam equation. A flexible beam element will be described by at least four 
partial differential equations. 


Lateral Bending 


The beam equations represent (1) Euler bending stiffness, (2) axial force stiffness, 
and (3) Torsion. For bending in the x-z plane: 

+ EI xJ u xXzzzz +GAu xXz2tt +F 0 u xXzz + K x (u xl + u xX ) = F xd {z,t) 

F x,iM ( 10 ) 


Axial Deformation 

Axial dynamics is represented by a wave equation with an additional term which 
represents a spring connected to a second distributed mass. 

mu zd - EAu zz l + K z (u z 2 + u zd ) = F zd (z,r) 


Torsion 

Torsional dynamics is represented by a wave equation 

ly, l^V,l — ^y,l w ycz,l ■*" Ky(lly/'2 + u y/,l) = ^z, l( z »0 ^2) 

Solution of the Partial Differential Equations 

The solutions of these partial differential equations for zero damping produce the 
sinusoidal and hyperbolic spatial equations which comprise the mode shape functions. For 
the case that F 0 = 0, the bending mode shape in the x-z plane is: 
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u x (z) = A x sin P lx z + B x cos p lx z + C x sinh/3 ljt z + D x cosh (5 lx z 


(13) 


Similarly, for bending in the y-z plane: 

u y (z) = A y sin fiiyZ + By cos p ly z + C y sinh /5 ly z + D y cosh $\ y z 


(14) 


The undamped mode shape functions for torsion and elongation about the z axis 
arc: 


u z (z) = A z sin/3 z z + B z cos/5 z z 
Uy(z) = Ay SmfiyZ + By COS/5 yZ 


(15) 

(16) 


These undamped mode shapes are expected to be good approximations to the exact 
solutions for low level of damping. The mode shape of the entire configuration consists of 
these functions, repeated for each beam element Because bending in two directions, 
torsion and elongation are considered, a total of 12 coefficients arc needed. The vector of 
coefficients is the state vector of the structural dynamics. A vector of the coefficients of 
these sinusoidal and hyperbolic functions will serve as the state vector. 


0 T = [a x B x C x D x Ay By Cy Dy A z B z Ay By] 


(17) 


Under conditions of applied forces it is necessary to include rigid body modes. 
Their coefficients will expand the state vector accordingly. All deflections, forces, 
moments, and accelerations will be expressed in terms of such state vectors. 

The motion of each rigid body is put in terms of the deflection at the point of 
attachment of a particular reference beam element The linear and angular deflection vectors 
can be expressed as: 


u = Q u (z)G 
u = Q s (z)Q 


(18) 

(19) 
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Next, it is necessary to express the forces and moments at either end of the beam 
elements. The force and moment vectors are: 


^attach - (20) 

^attach =/ m( z )© (21) 


It is also necessary to account for changes in axes from each beam to the body to 
which it is attached, and for points of attachment at some distance away from the center of 
gravity. The force and moment that a beam-i applies to a body-j are: 


^body-y ^body — j ^beam-i^V,/ (“)© 

^body- y — ^body- y ^beam- i(^b earn- ,•/>(*)}© 


( 22 ) 

( 23 ) 


The partial differential equations provide the relationships between the modal 
frequency and the eigenvalues for the mode shape equations. The lateral beam, axial 
deformation and torsion equations can be solved for the zero damping cases to produce the 
following relationships between the modal frequency and the wave numbers in the mode 
shape function. 


For bending in the x-z plane: 

p w =- 5b + Vcs^+wy 2 /*/, 

Au 2 = ~ 5b + ^(-5b) 2 +mco 2 / EI X 
where b = mco 2 / GA + F 0 / EI X . 


( 24 ) 

( 25 ) 


The case for bending in the y-z plane is similar. For torsion and elongation: 
Py = CO / ^Gly / m 

P 2 = co / ^EA~fm 


( 26 ) 

( 27 ) 
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Transfer Matrix Approach 


The transfer matrix approach [14] is suitable for large systems made up of several 
subsystems. The typical subsystem may be simple elements like a scalar spring or a 
complex Bemoulli-Euler or Timoshenko beam element The subsystems are cast in the 
form of a field and a point matrix. The formulation is in terms of the state vector which is a 
column matrix of displacements and internal forces. The treatment of the transfer matrix 
derivation for rigid bodies and flexible beams follows the work in [15]. 


Rigid Body 


The translational and rotational equations of the/^ body can be described by the 
following equations (figure 2). 


A 



Figure 2: Free-Body Diagram of the rigid body 

m jycM = Qf ~ Qf 

/ jf = Mj - Qfrf - Qfrf 


(28) 

(29) 


For harmonic motion, the equations are rewritten as 

Qf -Qf- m j^CM 

Mf = Mf + Qfrf + Qfrf - IjC0 2 Yj L 


(30) 

(31) 


The displacement of the center of mass Ycm is related to Yf and Yf by 
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(32) 

(33) 


1 f-l'ew-'fiy* 
rf=l cu + rfyj L 


Using the expression for Ycm from equation (33) and by the property of slope 
continuity, equation (32) can be rewritten as 


yi? _ yZ/ y/L 

*j ” l j r i l i 


(34) 


where r } = rf + rf . 

Substituting for *cm in equation (30), we get 

Qf=Qf-m J o> 2 (Yf--r i r; L ) 

£ 

Substituting the expression for Qj in equation (3 1), we get 

M? = Mf + rjQ f - mjdp-rfYf - (lj - m jr frf)(0 2 Yj L 


(35) 


(36) 


Equations (34), (35), (36) and the slope continuity condition yield the point matrix 
[PM]j for the rigid body element 


[PM]j = 


1 

0 


-ntjO) 

,.2 

rri j CO Vj 


-r } 0 O’ 

1 0 0 

mjCd 2 rj' 1 0 

-(ij-mfjbf) ij 1 


(37) 


’Y~ 

R 

~Y~ 

Y’ 

ii 

Y' 

Q 

Q 

M 

j 

M 


( 38 ) 
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Flexible Beam 


The field matrix for Bemoulli-Euler beam is derived from the solution of the 
bending mode slope (Equation 13) as follows. 


At the left end of the beam (z = 0) , the displacement Y (0) slope T'(0) , shear force 
Q(0) and bending moment M (0) will be 


~Y~ 

R 

■ 0 

1 

0 

1 ■ 

"V 

r 


p 

0 

p 

0 

B x 

Q 


-kp 3 

0 

+kp 3 

0 

c x 

_M_ 

h 1 

0 

-kp 2 

0 

-kP\ 

Px. 


(39) 


where Q = kY"' and M = kY",k = EL For notational simplicity, the subscript* 
on P is dropped. At the right end of the beam (z = L) 


~Y~ 

L 

sin PL 

cos PL 

sinh PL 

coshpL 

X' 

Y' 


P cos PL 

-p sin pL 

P coshpL 

PsinhpL 

B x 

Q 


kp 3 cos pL 

-kp 3 sir \pL 

kp 3 cosh PL 

-kp 3 sinh pL 

C x 

M 

j 

_kp 2 sin pL 

kp 2 cos PL 

-kp 2 sinh pL 

kp 2 coshpL _ 

Px. 


Solving for the coefficients A X ,B X ,C X ,D X from equation (39) and substituting in 
equation (40) we get the field matrix for the beam element as: 


[FM] r 


■^•(cos/JL + cosh/JL) 
-Xsin pL + sinh fiL) 
— (sin PL + sinh pL) 
^-(cos pL + cosh/JL) 


— (sin PL + sinh pL) 
^•(cos/3L+cosh/JL) 
^kp 2 (cosPL + cosh/JL) 
^■(sin/JL + sinh/iL) 


^^{smpL-smhpL) 

2 ^ (cos pt - cosh pL) 

y- (cos PL - COSh.pL) 

■^-{sinpL-swhpL) 

ip 


( cos PL ” cosh P L ) 
j^{smpL-swhpL) 
PL- sinh pL) 
^•(cos/JL -coshpL) 


’Y~ 

L 

~Y~ 

Y' 

1 

II 

Y' 

Q 

Q 


j 

M 


(41) 
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The transfer matrix for a flexible beam with a mass at the right end is 



Figure 3: Beam Offset 
Offset Attachment 

The planar offset attachment transfer matrix can derived from figure 3. The offset 
of point 2 from the origin (point 1) is given by r x and r z . 


X' 


"i 

0 

~ r z 

0 

0 

O' 

X" 

« 2 


0 

1 

r x 

0 

0 

0 

u z 



0 

0 

1 

0 

0 

0 

°y 

F x 


0 

0 

0 

1 

0 

0 

F x 

Fz 


0 

0 

0 

0 

1 

0 

F z 

My 

2 

0 

0 

0 

r z 

~ r x 

1_ 

My 
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Joint with Compliance 


The joint compliance transfer matrix is derived for planar motion assuming a spring 
of stiffness k Q . The joint equations for planar motion are 


K) 2 = K), 

[M,\=k e (e i-Oj) 


(44) 

(45) 


The transfer matrix can be expressed as 




’i 

0 

0 

0 

0 

0 ' 


«Z 


0 

1 

0 

0 

0 

0 

«z 



0 

0 

1 

0 

0 

-m e 


Fx 


0 

0 

0 

1 

0 

0 

P X 

Pz 


0 

0 

0 

0 

1 

0 

Pz 

My 

2 

_0 

0 

0 

0 

0 

1 

< M yJ i 


(46) 


Rigid Body Control 


The rigid body point matrix for a body with mass and inertia but with r ,• equal to 
zero is obtained from equation (37). Using the Laplace variable s in the place of -to , 
we get 




■ 1 

0 

0 

0 

0" 




0 

1 

0 

0 

0 

«z 



0 

0 

1 

0 

0 


Fx 


Ms 2 

0 

0 

1 

0 

Fx 

Pz 


0 

Ms 2 

0 

0 

0 

Pz 

w 

2 

_ 0 

0 

Is 2 

0 

1_ 

My 

yj i 


(47) 


Rigid body controllers basically stabilize the system asymptotically and are of the 
proportional derivative type. The transfer matrix is modified in the following manner for 
control: 
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Joint Control 


ms 2 — » ms 2 + c x s + k x 
Is 2 -» Is 2 + c 9 s + k e 


Analogous to the rigid body controller, it is also possible to embed local feedback 
control effects into the transfer matrix for a joint Only one matrix element is affected. For 
proportional derivative control including sensor f(s) and actuator g(s) dynamics, the matrix 
in equation (46) is modified thus: 

Hkg^l/ [/(j)g(jX*t + M] ( 48 ) 


Alignment Matrix 


The alignment of each element with respect to global coordinates is accomplished 
by a simple matrix multiplication. The planar alignment transfer matrix is 


x~ 


cos a 

-sin a 

0 

0 

0 

O' 

X' 

«z 


sin a 

cos a 

0 

0 

0 

0 

“z 

*7 


0 

0 

1 

0 

0 

0 


F, 


0 

0 

0 

cos a 

-sin a 

0 

F x 

Fz 


0 

0 

0 

sin a 

cos a 

0 

F z 

m 2 _ 

2 

0 

0 

0 

0 

0 

1 

M z 


End-to-End Transfer Matrix 


(49) 


The transfer matrix which relates the deflections and loads at the space shuttle to 
those through the RMS to the Space Station consists of the product of all of the elements as 
shown in Figure 4. 
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Figure 4: Shuttle/RMS/Station Configuration 


'Y‘ 




~Y~ 

Y' 

Q 


[^station ^wrist ^boom 2 ^elhow 

^boom 1 ^shoulder ^shuttle] 

Y' 

Q 

M 

station 



M 


J shuttle 


•u t 


total 


Y' 

r 
Q 

L^Jshuttle 


(50) 


It is now possible to derive the characteristic equation for the total system. The 
effect of all the control systems will be reflected in the characteristic equation since they 
form a part of the rigid body and joint transfer matrices. 

The transfer matrix forms an intermediate step in the computation of the 
characteristic equation. For beams and masses connected to one another, the transfer 
matrix between station 1 and station n is derived by multiplying the appropriate field and 
point matrices. The expression for a typical problem may be expressed as 


~Y~ 


011 012 013 014 

■y 

Y' 


021 022 023 024 

Y' 

Q 


031 032 033 034 

Q 

M 

n 

1 

•e- 

m 

(N 

Tf 

1 
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If the boundary is cantilevered, then 


m Y~ 


‘O' 

r 

= <D 

0 

0 


Q 

_o_ 

n 

_M 


( 52 ) 


Rearranging the state vector we have 


[A] 


X 

y 1 

L n 

Q i 

l Mi. 


= [ 0 ] 


(53) 


where 


[A] = 


'-1 

0 

013 

014 

0 

-1 

023 

024 

0 

0 

033 

034 

0 

0 

043 

044 


(54) 


The characteristic equation is given by 

det[A] = 0 ( 55 ) 

or 

det[<l>22] = 0 


or 

033 044 ~ 034 043 = ® (56) 

For continuum models, equation (56) has infinite solutions and is solved by search 
techniques to determine the frequencies. Similar characteristic equations can be derived for 
other boundary conditions. 
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For free-free boundary conditions, the characteristic equation is given by 


det[d>2i] = 0 


The transfer function relating the linear and angular deflection of the shuttle (for 
example) to applied forces and moments are: 


~Y~ 

Y' 


— <J>2i <$22 


Ishuttle 


F 

M 


Jshuttle 


The Shuttle/RMS/Payload configuration is now studied from the continuum 
viewpoint. 


Study of Shuttle/RMS/Payload Assembly 

The Space Shuttle/RMS/Payload assembly is modeled and analyzed using the 
continuum and the finite element approach. For the continuum analysis, the planar transfer 
matrix approach is used to generate the frequencies of the configuration and the transient 
response of the structure. 

The data for the two link RMS configuration is extracted from the payload 
deployment and retrieval document [16]. Links 3 and 4 of the RMS arm are used in the 
simulation. In this work, each link was assumed to be made-up of one material with 
uniform section properties unlike reference [16] where the links were made-up of 3 
segments each with different properties. The link properties are listed in Table 1 . 

The space shutde and payload are modeled as rigid bodies with a mass of 6176 
slugs and 124.22 slugs respectively. The inertia I yy of the space shuttle is 6.99 E6 lbs-in . 





Payload 



Figure 5: Shuttle/RMS/Payload Configuration 

Figure 5 shows the Shuttle/RMS/Payload Configuration. The following three cases are 
considered in this work: 

Case 1: Shuttle with zero inertia and offset 

Case 2: Shutde with inertia and zero offset 

Case 3: Shuttle with inertia and offset 

For the transfer matrix approach, the relationship between the shuttle and the 
payload is 

[7iV/] — |t^p a yi oa( j Cjjnk 2 ^]j n k i ^shuttle] 

The characteristic equation for the free-free configuration is derived and the 
frequencies are evaluated. In order to obtain the y-z bending frequencies, the characteristic 
equation is again solved numerically using the appropriate flexural rigidity value. The 
results are compared with the frequencies from PDEMOD. 

The NASTRAN model of the Shuttle/RMS/Payload assembly consisted of the RMS 
being modeled using 50 bar elements each. The shuttle and the payload were modeled 
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using two concentrated masses at either end of the RMS. For the shuttle with inertia and 
offset case, the concentrated mass card in NASTRAN was suitably modified. 

Table 1: RMS Link Properties 


Property 

Link 1 

Link 2 

mass 

9.5485 slugs 

5.9901 slugs 

length 

21.0 ft 

23.0 ft 

Eiyy 

5.6458 E6 lbft 2 

3.4166 E6 lbft 2 

F.T7.7. 

5.2083 E6 lbft 2 

2.4375 E6 lbft 2 

G 

3.846 E5 psf 

3.846 E5 psf 

J 

6.7711 ft 4 

5.0558 ft 4 

m/1 

0.4547 slug/ft 

0.2604 slug/ft 

A 

0.9218 ft 2 

0.9218 ft 2 


The frequency spectrum in Hertz of the three configurations is shown in Tables 2-4. 


Table 2: Frequencies for Case 1 


Mode 


NASTRAN 


PDEMOD 


Transfer 

7 

2.528 

2.862 

10.629 

11.747 

? 

23. ni 

25.9641 


1 

2 

3 

4 

5 

6 

7 

8 


2.054 

2.527 

2.862 

10.629 

11.747 

18.479 

23.171 

25.963 


2.528 

2.863 

10.63 

11.74 

? 

23.17 

25.96 





Amplitude 


Table 3: Frequencies for Case 2 


Mode 

NASTRAN 

PDEMOD 

Transfer 

1 

0.1825 

0.1961 

” 0.1892 

2 

2.0548 

2.0133 

? 

3 

2.8627 

? 

? 

4 

4.2036 

4.2065 

4.2033 

5 

? 

4.6983 

4.6973 

6 


7.6745 

? 

7 

11.7471 

? 

? 

8 

13.3324 

13.3356 

13.3337 


Table 4: Frequencies for Case 3 


Mode 

NASTRAN 

Transfer 

1 

? 

0.5428 

2 

? 

0.3685 

3 

4.2037 

4.4246 

4 

4.7001 

4.9433 

5 

? 

9.8986 

6 

? 

9.9017 

7 

13.3339 

13.5996 

8 

14.7421 

15.0281 


Figure 6 shows a transient response obtained from PDEMOD for a similar configuration 
with and without joint control. The results show the promise of the continuum approach. 



Figure 6: Transient response of MB-1 Configuration 
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CONCLUDING REMARKS 


Partial differential equation models of flexible structures offer significant 
advantages over finite element models for parameter estimation and control studies because 
of the smaller number of model parameters. Until recently work was needed to generate 
distributed parameter models of complex configurations which were also flexible. The 
computer program, PDEMOD, enables the generation of distributed parameter models of 
flexible spacecraft Any configuration which can be modeled by a network of flexible 
beam elements and rigid bodies can be modeled using PDEMOD. The modeling process is 
well suited for the evolving Space Station Freedom, for the cases in which (1) the Space 
Station assembly is attached to the Shuttle, (2) the assembly is linked to the Shuttle through 
the RMS arm, and (3) the Space Station assembly is free of the Shuttle. 

Comparisons of the model accuracy of finite element and continuum models of 
flexible structures point out the limitations of finite element modeling. First, the level of 
complexity that is practical for finite element models is limited because of the computational 
burden. The result is a limit to the accuracy that can be obtained. Second, as high levels of 
accuracy are sought using finite element models, the difficulties in solving the eigenvalue 
problem become more significant It is quite possible, then, that for certain applications 
continuum models can be more accurate. 

A distributed parameter model of the Space Shuttle-RMS was generated using the 
transfer matrix method and the software PDEMOD. The results show a very good 
agreement with a detailed finite element model. Future directions include the frequency 
characterization of structures with embedded control. 
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